arXiv: 1502.02979v2 [astro-ph.HE] 1 Apr 2015 


Magnetar superconductivity versus magnetism: Neutrino cooling processes 

Monika Sinha^’0 and Armen Sedrakian^ 

^Institute for Theoretical Physics, J. W. Goethe-University, D-60A38 Frankfurt am Main, Germany 

(Dated: April 2, 2015) 

We describe the microphysics, phenomenology, and astrophysical implication of a _B-field induced 
unpairing effect that may occur in magnetars, if the local B-held in the core of a magnetar exceeds 
a critical value Hc 2 - Using the Ginzburg-Landau theory of superconductivity, we derive the Hc 2 
held for proton condensate taking into the correction (< 30%) which arises from its coupling to the 
background neutron condensate. The density dependence of pairing of proton condensate implies 
that Hc 2 is maximal at the crust-core interface and decreases towards the center of the star. As a 
consequence, magnetar cores with homogenous constant helds will be partially superconducting for 
“medium-held” magnetars (10^® < B < 5x 10^® G) whereas “strong-held” magnetars (B > 5x 10^® 

G) will be void of superconductivity. The neutrino emissivity of a magnetar’s core changes in a 
twofold manner: (i) the B-held assisted direct Urea process is enhanced by orders of magnitude, 
because of the unpairing effect in regions where B > Hc 2 ', (ii) the Cooper-pair breaking processes on 
protons vanish in these regions and the overall emissivity by the pair-breaking processes is reduced 
by a factor of only a few. 


I. INTRODUCTION 


The protonic fluid in the cores of magnetized neutron 
stars is a type-II superconductor, i.e., it supports the 
magnetic B field by forming quantized electromagnetic 
vortices with density n„ = B/^q, where $o = 7rhc/e 
is the quantum of flux, in the field range Hd < B < 
Hc 2 [H- The lower critical field Hd is the held strength 
at which the emergence of the hrst vortex (hux tube) 
becomes energetically favorable. At the upper critical 
held strength Hc 2 the normal cores of the vortices touch 
each other and superconductivity is destroyed. The Hc 2 
held is density-dependent and is given by 


H,2 = 


^>0 

2 <’ 


( 1 ) 


where is the coherence length of the proton conden¬ 
sate, which scales inversely with the pairing gap A. The 
coherence length appears in Eq. (HD because Hc 2 is the 
held at which the Larmor radius of protons becomes com¬ 
parable to the size of a Cooper pair ~ . A held B ^ Hc 2 
disrupts the coherence among the protons which form a 
Cooper pair and, therefore, destroys their superconduc¬ 
tivity. For the proton superconductor in the cores of neu¬ 
tron stars 10^® < Hc 2 < 10^^ G, i.e., Hc 2 is well above 
the helds expected in the interiors of ordinary neutron 
stars {B ~ q). 

The inferred magnetic helds on the surfaces of magne¬ 
tars are of the order of 10^® G. Their interior B helds 
are not known, but it has been frequently conjectured 
that they are larger than the surface held. The conjec¬ 
tured maximal B held, which is consistent with the virial 
theorem for self-gravitating magnetic equilibria, is of the 


* Present address: Indian Institute of Technology Rajasthan, Old 
Residency Road, Ratanada, Jodhpur 342011, India 


order Bmax — 10^® G. Because i?max > Hc 2 and be¬ 
cause these helds may vary over the star’s core, we may 
anticipate an intimate interplay between the magnetism 
and superconductivity in the interiors of magnetars de¬ 
pending on whether the local held is above or below Hc 2 - 
Some observational arguments where put forward in re¬ 
cent years in favor of type-I superconductivity i. 
However, our choice of the equations of state of dense 
matter and microscopic parameters of the proton super¬ 
conductor predict type-II superconductivity throughout 
most of the core of a neutron star, as we show below. 

The purpose of this work is to show that large enough 
magnetic helds in the interiors of magnetars unpair pro¬ 
ton superconductor in a strongly density-dependent man¬ 
ner. We then go on to study the consequences of this 
magnetically induced unpairing effect on the neutrino 
emissivity of magnetars. Neutrino emissivities are the 
key ingredients for the simulations of thermal evolution 
of magnetars and can be confronted with the measured 
x-ray huxes from the surfaces of magnetars. 

This paper is structured as follows. In Sec. Ullwe dis¬ 
cuss the input physics, i.e., the underlying equation of 
state (EoS) and composition of matter which sets the 
stage for the following discussion. In Sec. uni starting 
from the Ginzburg-Landau (GL) functional for proton 
superconductor coupled to neutron superhuid, we derive 
an expression Hc 2 , which accounts for the density-density 
coupling between the proton and neutron condensates. 
In Sec. IIVI we compute the neutrino emissivities of the 
Urea process and pair-breaking processes in magnetars 
including the unpairing effect. Our conclusions and an 
outlook are given in Sec. |Vl 


II. MICROPHYSICAL INPUT 

Consider a magnetar with a nonstrange baryonic core 
consisting of neutrons (n), protons (p), electrons (e), 
and muons (p) in /3 equilibrium. We choose to work 
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with a relativistic density functional (DF) with density- 
dependent couplings derived in Ref. 0 to obtain the 
equation of state (EoS) and composition of matter in the 
star’s core and inner crust. The latter parametrization is 
in excellent agreement with the nuclear phenomenology 
as it predicts saturation density no = 0.152 fm“^, bind¬ 
ing energy per nucleon E/A = —16.14 MeV, incompress¬ 
ibility Kq = 250.90 MeV, symmetry energy J = 32.30 
MeV, symmetry energy slope L = 51.24 MeV, and sym¬ 
metry incompressibility Ksym = —87.19 MeV all taken 
at saturation density [3- For completeness we show the 
EoS of baryonic matter in Fig. [1] derived from this DF. 
Compact star models based on this DF were constructed 
elsewhere @ where it has been shown that the resulting 
maximum mass predicted by this EoS is well above the 
current observational lower limit 2Mq on the maximum 
mass of any compact star. Strangeness in form of hyper- 
ons or deconfined two- or three-fla.vor mia.rk matter mav 
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FIG. 1. Zero-temperature equation of state of dense mat¬ 
ter composed of neutrons, protons, electrons, and muons 
in P equilibrium derived from the relativistic DF with the 
parametrization of Ref. [j. 


The composition of dense matter corresponding to our 
EoS is shown in Fig. [21 where we show the abundances of 
species n^/nb, where i £ n,p, e, /i as a function of baryon 
density rib normalized by the nuclear saturation density 
of the DF. The abundances of protons and electrons are 
equal up to the point where muons set in. The threshold 
value of Urea process in non-magnetized matter np/nt ~ 
0.11 is reached at the density n ~ 3no. The composition 
of matter itself will be affected by a strong B field, when 
electromagnetic interactions become of the order of the 
nuclear scale set by the Fermi energies of the constituents. 
However, below the field values 10^® G the abundances of 





FIG. 2. (Golor online) Dependence of particle abundances 
Ui/rib, i € n,p,e,fj, on the net baryon density nj in units of 
saturation density no — 0.152 fm“®. The vertical line shows 
the approximate Urea threshold Vurca = 0.11 for proton frac¬ 
tion. 


baryons for non-zero B are indistinguishable from those 
in the B = 0 case (see Ref. and references therein). 

The pairing channels in neutron star matter corre¬ 
spond to the attractive most dominant phase-shifts at 
given density or energy of nucleons (for a review see, 
e.g.. Ref. Q). Low density neutron matter in the crust 
of compact stars pairs in the channel; above the sat¬ 
uration density the neutron fraction is large enough (and 
energies are high enough) to render the S'-wave interac¬ 
tion repulsive. The dominant pairing at these densities is 
in the channel. Protons are in the continuum in 

the fluid core of the star and are much less abundant than 
neutrons, therefore their energies are low enough to fa¬ 
vor the pairing. The pairing gaps in these dominant 
channels adopted from Refs. are shown in Fig.jSl The 
formulas which Ht these gaps are listed in the Appendix. 
At asymptotically high densities pairing of protons in the 
^P 2 -^F 2 channel may occur. In the exceptional models 
where matter is nearly isospin symmetrical at high den¬ 
sities, spin-one, isospin-zero pair ing in the ®Zl 2 channel 
may become the dominant one [lOj. Examples of such 
models are those which feature kaon condensates. Below 
we do not consider proton P-wave pairing or neutron- 
proton D-wave pairing. 


III. GINZBURG-LANDAU THEORY OF 77^2 IN 
DENSE MATTER 

Type-II superconductivity is characterized by the GL 
parameter, k = where Sl is the London penetra- 
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FIG. 3. (Color online) Dependence of S- and P-wave pairing 
gaps of neutrons (dash-dotted and solid lines) and of S-wave 
protons (dashed line) on their respective Fermi momenta. 


tion depth of the B held in a superconductor, having the 
range 

-^ < K < oo. (2) 


The critical value Kc = = 0.7071 separates the 

domains of type-I and type-II superconductivity. 

The magnetic held is conhned to electromagnetic vor¬ 
tices for held values between the lower Hd and upper Hc 2 
critical helds. If the B held is lager than Hc 2 it unpairs 
the Cooper pairs and the material makes a transition to 
the normal state. The phase transition from supercon¬ 
ducting to the normal state in the vicinity of Hc 2 can be 
described in terms of the GL theory, because the super¬ 
conducting order parameter is small. Note that in the 
vicinity of Hc 2 the superconducting order parameter is 
small because of the large B held and the GL expansion 
is valid not only near the critical temperature Tc, but in 
the entire temperature range 0 < T < T^. 

We start by writing down the GL functional for a su- 
perhuid neutron and superconducting proton mixture 


F[(j),tjj] = P„[(/)] -P aT\tjjf + 


1 

4mp 



Stt ’ 


( 3 ) 


where ijj and (j) are the proton and neutron condensate 
wave-functions, nip is the proton mass, t = {T—Tcp)/Tcp, 
where T^p is the critical temperature of superconducting 
phase transition of protons. Here a and b are the coeffi¬ 
cients of the GL expansion for the proton condensate, b' 


describes the density-density coupling between the neu¬ 
tron and proton condensates. This type of GL functional 
was analyzed initially to study the current-current cou¬ 
pling between the neutron and proton condensates m 
(the entrainment effect, see Ref. M)- More recent 
study of Ref. Q discusses the density-density coupling 
between the neutron and proton condensates and pro¬ 
vides the relevant microscopic expressions for the coeffi¬ 
cient of the GL functional. The effective vector potential 
can be decomposed as A = Agm + ^ent, where the first 
term is the ordinary vector potential of electromagnetism 
with B — 'V X A. The second term is the “entrainment” 
vector potential Agnt = (/tc/e)[(m* —mn)/mn] V((), where 
TO* is the proton effective mass, see Ref. [Il|. The en¬ 
trainment effect describes the current-current coupling 
between the neutron and proton condensates. 

The explicit form of the contribution of the neutron 
condensate to the GL functional, Fn[<p], is not required 
in the following. 

The minimization of the GL functional with respect to 
tp*, i. e., 6F[(j),'il!]/6ip* = 0 gives 

— —A^ Ip + arip + b\p)\^ip -|- b'\(p\^p) = 0. 

( 4 ) 

The equilibrium value of the condensate is given by the 
solution of Eq. O 

ipi^ar+ b\'tp\'^+ b'\(p\^) = Q, (5) 

from which we obtain the two possible equilibrium solu¬ 
tions 


= 0, T > Tg, (6) 

\P^\^ = -\{aT + b'\cP\% T<n. (7) 

The variation of the GL functional with respect to the 
electromagnetic vector potential 6F[(p,tp]/6A = 0 gives 

±VxVxA = j, (8) 

where 

3 = - {lp*Vlp - IpVlp*) -IV'I (Aem + 2Aent) 

m me 

(9) 

is the proton super-current. It contains the conventional 
electromagnetic current oc V'i/; as well as the entrainment 
current oc Aent oc V(/). 

Equations (|1]), ([5]) , and ([HI) constitute the GL equations 
in their most general form. To derive the value of the 
upper critical field i7c2 it is sufficient to keep only the 
linear in tp terms in the GL equations above. To this 
order Eq. reduces to 

V xV X A = 0 + O{\tp\^). (10) 

Furthermore, because V x = 0 identically (except 
at the singular points where the neutron vortices are 







located), we can make the replacement A —>■ Aem in 
Eq. (llOp . The small-scale (local) magnetic field is ho¬ 
mogenous, therefore the corresponding vector potential 
Aem is linear in coordinates. We choose Aom along one 
of the directions of the Cartesian system of coordinates, 
say z direction, without loss of generality. Assume that 
the B field is in the y direction. Then, ip = ip{x) only. 
To linear order in ip the solution of Eq. is Aem = Bx. 
Substituting this into the first GL equation (jd]) one finds 

“ + o(IV'P)- 

( 11 ) 


The mathematical form of this equation is that of the 
harmonic oscillator, therefore, its solutions is read off as 



4,ttB 

$0 


( 12 ) 


We are interested in the strongest field for which solutions 
with Ip ^ 0 are still possible. This is the case n = 0 
in Eq. m which identifies the critical field B = Hc 2 - 
Consequently 


Hc2 


27r 


$0 


ArUp 


(ar -I- b'\(pf) 


1-H 


\b'm^ 


ar 


(13) 


where we used the relation (mplarj)^/^ = h/2^p and the 
fact that 6' < 0, see below. If 6' = 0 Eq. (US reduces 
to the standard result [I^ . To evaluate the correction to 
the critical field note that | q ; t | = \ipo\'^b and, therefore, 

{b'\p\'^/aT) = {nn/np){\b'\/\b\). 

The coefficient b' which takes into account beyond 
mean-field coupling between the neutron and proton 
condensates was computed by Alford et al. Q in /3- 
equilibrated, charge-neutral nuclear matter diagrammat- 
ically. They also provide the mean-field expression for b. 
Using their results we find that 


n„ |b'| _ 277r^ nl Ag 
Up \b\ 4 ylulmpkF^' 


(14) 


where we used the value of b' valid in the regime Ap <g; yp 
and A„ <g; and T —?► 0 and the value of parameter 
Gnp = 10“® MeV“^ 0. Note that close to the critical 
temperature Tc alternative expressions provided by Al¬ 
ford et al. 3 should be used. The correction in Eq. m 
owing to the coupling between the neutron and proton 
condensates is < 30%; it is small because the coupling be¬ 
tween the condensates arises only via fluctuations which 
vanish in the ground state. The main uncertainty in 
Eq. (fHl) is the contact pairing interaction in the isosinglet 
channel Gnp', its value quoted above should be viewed as 
an order of magnitude estimate. An additional uncer¬ 
tainty arises from the not well-known value of the gap in 
the proton spectrum Ap which may vary by a factor of 
few. 



FIG. 4. (Color online) Upper panel: Dependence of pairing 
gaps for nentrons (^So and ^P 2 channels) and for protons {^So 
channel) on baryonic density in units of nuclear saturation 
density. Lower panel: Dependence of the critical unpairing 
field Hc 2 on baryonic density with account for the coupling 
between the neutron and proton condensates (full line) and 
without (dashed line). 


The analogy between Eq. m and the one describ¬ 
ing harmonic oscillator in quantum mechanics can be ex¬ 
ploited further to write down the most-general “harmonic 
oscillator” type solution of Eq. (EH), which describes a 
vortex in the x-y plane (with the field directed in the z 
direction). The corresponding wavefunction can be writ¬ 
ten as 


ip{x,y) = Gnexp[—KB(x—k/KB)^/2 + iky], (15) 

n— — co 

where the coefficient Gn and k depend on the type of the 
proton vortex lattice. Assuming triangular lattice one 
finds k = K(7r-\/3)^^^ and the set of condition C '„+4 = Cn, 
Go = Gi = G, G 2 = Co = —C, where C is given by 
the normalization of the wave-function to the density of 
condensate. 

Table I lists the key parameters of the proton supercon¬ 
ductor for a range of densities corresponding to the star’s 
fluid core. The coherence length has a minimum, which 
reflects the density dependence of the gap (^p oc A“^). 
Because Hc 2 scales inversely with ^p, the critical field has 
a maximum, with max Hc 2 = 7.37 x 10^® G at nt, = 0.7no 
in our setup. The London penetration depth scales as the 
root of inverse proton density, therefore it decreases as 
the density increases. This has the consequence that the 
GL parameter drops below the critical value and the 
proton superconductor becomes type-I. However, this oc- 
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rib/no 

kpp 

Ap 

m*/mp 

ip 

5l 

K 

Rc2 

0.140 

0.12 

0.02 

0.93 

76.1 

929.2 

12.2 

0.06 

0.300 

0.20 

0.24 

0.89 

11.9 

425.0 

35.6 

3.15 

0.500 

0.28 

0.55 

0.85 

8.0 

238.6 

29.8 

7.37 

0.700 

0.36 

0.76 

0.81 

7.8 

161.1 

20.6 

7.08 

0.900 

0.44 

0.85 

0.78 

8.7 

119.5 

13.7 

5.15 

1.100 

0.51 

0.86 

0.76 

10.4 

93.9 

9.1 

3.39 

1.300 

0.58 

0.81 

0.74 

13.0 

75.2 

5.8 

2.06 

1.500 

0.67 

0.73 

0.71 

17.0 

61.0 

3.6 

1.18 

1.700 

0.74 

0.62 

0.70 

22.8 

51.2 

2.2 

0.64 

1.900 

0.81 

0.45 

0.68 

35.0 

44.3 

1.3 

0.27 

2.100 

0.88 

0.16 

0.67 

106.4 

39.2 

0.4 

0.03 


TABLE I. Microscopic parameters of proton superconductor 
and the upper critical field for unpairing Hc 2 for a range of 
matter densities. 


curs only in the high-density end of the proton supercon¬ 
ductivity domain and should be relevant only for compact 
stars with central densities exceeding this value. 

FigureH] (upper panel) displays the dependence of pair¬ 
ing gaps on baryon density for the composition of matter 
implied by our chosen EoS. (Note that gaps displayed in 
Fig. Oas functions of Fermi momenta of particles are EoS 
independent, whereas those in Fig. |4] are specific to our 
EoS). The dependence of the Hc 2 held on density is shown 
in Fig. |4] (lower panel). It is seen that magnetars with 
interior fields with B < max Hc 2 — 7.37 x 10^® G will 
be partially superconducting, which means that regions 
where B < Hc 2 will be superconducting whereas the re¬ 
gions where B > Hc.2 are not. Clearly, magnetars with 
B > max Hc2 will be fully non-superconducting. The 
maximum of Hc 2 is attained close to the crust-core inter¬ 
face (corresponding to rib = O.Sno)- This implies that for 
partially superconducting magnetars with B < max Hc 2 
the unpairing by the magnetic field will remove proton 
superconductivity in the inner core, whereas the outer 
core could be still superconducting provided the B field 
is approximately homogeneous and constant in the fluid 
core of the star. This is a reasonable assumption, because 
the density gradients are small in the fluid core. 


IV. NEUTRINO EMISSIVITY OF MAGNETAR 
CORES 

This section studies the implications of the unpairing 
effect, discussed in the previous section, on the neutrino 
emission processes from dense matter in magnetars. We 
focus below on the neutrino emission processes which are 
dominant below the critical temperature Tcp of proton su¬ 
perconductivity, specifically the B field assisted Urea and 
the pair-breaking processes. The implications of the un¬ 


pairing effect for processes such as the modified Urea pro¬ 
cess and the modified bremsstrahlung process are analo¬ 
gous to those for the direct Urea process and the imple¬ 
mentations in numerical codes should be straightforward. 


A. Direct Urea process 

The direct Urea process is kinematically allowed only 
above the threshold Yurca = np/nt > 11% in ordinary 
low-field compact stars, because for low proton concen¬ 
trations the energy and momentum conservation cannot 
be fulfilled simultaneously [3. Strong B fields change 
the phase-space of baryons. As a consequence, the di¬ 
rect Urea process is allowed even below the threshold 
Yurca To characterize the kinematics of the Urea 

process in a i? field it is convenient to introduce the pa¬ 
rameter [l^ 


kpn - ikpe + kpp? 

U2 ^^Fp 

'^Pn 


(16) 


where Npp = k‘^p/2\e\B is the number of Landau levels 
populated by protons. Thus, for a; > 0 the Urea process is 
forbidden in the low-field limit, but can become operative 
in strong magnetic fields. 

If, under such conditions, the Urea process operates 
at a fraction of its strength, it can still be an important 
factor in cooling the star’s core, because other processes 
are by orders of magnitude weaker. For a; < 0 the Urea 
process is allowed and the role of the magnetic field is 
to induce “de Haas-van Alfven” type oscillations in the 
emissivity of this process as a function of B field. 

The second effect of the strong magnetic field on the 
Urea process (not discussed so far) is the effect of unpair¬ 
ing of the proton superconductor by the field. Proton and 
neutron pairings restrict the phase space available for the 
process and, as a consequence, the rate of the direct Urea 
process is suppressed. This suppression at asymptoti¬ 
cally low temperatures is given simply by an exponential 
quenching factor exp(—A/T) for each participating nu¬ 
cleon, where A is the relevant pairing gap, T is the tem¬ 
perature (more accurate treatments are given, e.g., in 
Ref. [13). As outlined in Sec. imi large B fields unpair 
the proton superconductor, therefore the suppression of 
the Urea neutrino emission due to the gap in proton 
quasiparticle spectrum will be absent, i.e., only neutron 
pairing will contribute to the suppression. Because the 
gap for neutrons in the P-wave channel is smaller than 
the one in the 5'-wave channel for protons (see Fig. 0]), 
the onset of suppression will strongly deviate from the 
one expected in the case of superconducting protons. 

We now illustrate these qualitative arguments by nu¬ 
merical examples. In our setup the proton fraction re¬ 
mains below Yurca in the density range where proton S- 
wave superconductivity exists, i.e., densities n < 3no, 
therefore we explore first the domain a; > 0, where Urea 
process is forbidden in the zero-field limit. 
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FIG. 5. (Color online) The emissivity of the Urea process in 
the forbidden region {x > 0) in units of the zero-field emissiv- 
ity €o at fixed density n = tiq for temperatures T = 0.01 and 
0.1 MeV. The Urea emissivity is shown for the cases of (a) 
normal matter (dotted line), (b) paired neutrons and normal 
protons (dashed lines) and (c) paired neutrons and protons 
(dashed-dotted lines). However the case (c) cannot be real¬ 
ized because of the unpairing effect for all the plotted val¬ 
ues of 5 > Hc 2 = 16.57. Note that for fixed density the 
scaling of the kinematical factor x along the B-axis is given 
through its dependence on the number of Landau levels, i.e., 
X oc oc 


The Urea emissivity for B 0 is written as [l^ 

(1 + (17) 

where Gp is the Fermi coupling constant, is the axial- 
vector coupling, are the effective masses of neutron 

and proton, /ie is the chemical potential of electrons, 
TZ function encodes modifications due to the field and 
Sn/p = exp(—A„/p/T) are the suppression factors arising 
owing to the pairing of neutrons (n) and protons (p). The 
quenching of proton superconductivity implies Bp = 1 in 
Eq. (1171) . To model the function TZ in the forbidden region 
we use an approximate polynomial fit to the functions 
shown in Fig. 1 of Ref. [^, which is given by logj^g TZ = 
-0.35942 - 0.506418X -b O.OlSOSOSx^ - 0.00140399x3. In 
the allowed domain we use the fit formula [l^ 

^ , cos 6 

“ 0.5816+ |x|i i92 ( ) 

where 0 = (1.211 -b 0.4823|x| -b 0.8453|x|2-333)/(i + 
1.438 |x| 3-^°®) which is valid in the range —20 < x < 0 
and Npp —>■ oo. Figure [S] displays the neutrino emissivity 
via the Urea process in the forbidden region as a func¬ 
tion of the B field at fixed density (n = no) and two 



logic B 


FIG. 6. The emissivity of the Urea process in the allowed 
region (x < 0) in units of the zero-field emissivity eo at fixed 
density n = 1.5 no for temperature T = 0.1 MeV. The mag¬ 
nitude of proton and neutron gaps are Ap = 0.73 MeV and 
An = 0.28 MeV. The emissivity is shown in the cases of (a) 
unpaired matter, (b) paired neutrons and for B > Hc 2 (un¬ 
paired protons) and (c) paired neutrons and for B < Hc 2 
(paired protons). 


values of temperature. The unpaired case coincides with 
the results of Refs. [H, [l^. Magnetic field allows the 
Urea process to operate with emissivity comparable with 
the emissivities of competing processes in the asymptot¬ 
ically large field region B —>• Bmax, as seen in Fig. O 
The pairing of neutrons and protons requires an addi¬ 
tional multiplicative factor SnSp = exp[—(A„ -b Ap)/r] 
in the neutrino emissivity. We show the cases Ap = 0 and 
Ap yb 0 assuming that the neutron pairing gap A„ yb 0 
and corresponds to its value at B = 0. Because for all 
B held values B > maxBc 2 (logic [maxBc 2 ] = 16.87) the 
unpairing effect requires Ap = 0; thus the case Ap yb 0 
is not realized physically, but provides a measure of the 
error of neglecting the unpairing effect. It is evident from 
Fig.[5]that this error is substantial and is the consequence 
of the fact that Ap ^ A„ in our example. This condition 
holds except at the edge of the density domain of interest, 
see Fig. [H Thus, the proton pairing, if allowed, would 
suppress the emissivity stronger than the neutron pair¬ 
ing, but because of the unpairing effect the Urea emissiv¬ 
ity is suppressed by the neutron superfluidity only. As a 
consequence the Urea emissivity would be enhanced from 
its value which neglects the unpairing effect. To explore 
the allowed region x < 0 of kinematics for the Urea pro¬ 
cess in strong magnetic held we have artihcially increased 
the Fermi momenta of protons and electrons by factor 
of two. (In our models the proton fraction exceeds the 
Urea threshold at density which is larger than the max- 
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imal density at which proton S'-wave superconductivity 
exists). We also choose to work at density 1.5 uq because 
Urea process becomes operative in the high density do¬ 
main. Figure [5] displays the neutrino emissivity of the 
Urea process in the allowed region as a function of the B 
field for unpaired matter and for cases B < Hc 2 (super¬ 
conducting protons) and B > Hc 2 (non-superconducting 
protons). In the case of unpaired neutrons and protons 
(panel a), the B field induces de-Haas-van Alfven type 
oscillations in the emissivit y ar ound its value in the zero 
B field limit, as expected [isl H^ . For fields B < 77^2 
the emissivity is suppressed by neutron and proton pair¬ 
ing simultaneously; for B > Hc 2 protons are unpaired 
and the suppression is only due to paired neutrons. The 
transition from one regime to the other is seen as a jump 
in the emissivity in (b) and (c) of Fig. [5] at 77 = Hc 2 - 
The oscillations in (b) are around a value of emissivity 
which is about an order of magnitude smaller than the 
emissivity in the normal state, which reflects the suppres¬ 
sion via neutron pairing only. If the unpairing effect was 
neglected the emissivity would have remained about 4-5 
orders of magnitude below B = 0 case. 

Strong magnetic fields will influence the neutron su¬ 
perfluidity in the curst (S'-wave) and in the core (P- 
wave) differently. The S-wave condensate forms spin-zero 
Cooper pairs and the Pauli paramagnetic alignment of 
neutron spins along the B field will act to quench their 
pairing. Generally, this quenching is effective for fields 
larger than those discussed here (P > 10^”^ G), but the 
value of the critical field depends on the gap in the zero 
field limit, which has an uncertainty of an order of mag¬ 
nitude. The P-wave condensate forms spin-one Gooper 
pairs and the magnetic field will align the spins of Gooper 
pairs without affecting their internal structure. Initial 
studies of P-wave pairing in strong fields show that there 
is no suppression of the pairing induced by the field [l^ . 
Therefore, as far as the Urea process is concerned, we do 
not expect additional suppression of pairing due to the 
B field in the P-wave paired core. 

Thus we conclude that the unpairing effect which de¬ 
stroys the proton condensate can strongly influence the 
neutrino emissivity via the Urea process in the cores 
of magnetars. These modifications may have important 
consequences on the modeling of thermal transients and 
cooling in magnetars. 

There exists an additional channel of neutrino losses, 
which arises once the interaction energy of the B field 
with the spin of a nucleon becomes of the order of temper¬ 
ature - the direct bremsstrahlung process N ^ N + 1 / + ^, 
where N refers to a nucleon. This process is strictly for¬ 
bidden in the non-magnetic case, but becomes operative 
in a strong enough B field, because of the paramagnetic 
splitting of the energies of nucleons in a strong B field. 
Spin-flip neutrino emission is effective within the win¬ 
dow of splitting energies of the order of the temperature 
of ambient matter [l9|. By the same argument as in the 
case of the Urea process above, the bremsstrahlung pro¬ 
cess p ^ p + v + D will remain intact in magnetars, in 


contrast to the case of ordinary neutron stars, where it 
would be suppressed by the proton pairing at low enough 
temperatures. 


B. Pair-breaking processes 


The formation of nucleonic BGS condensates leads to 
the pair-breaking neutrino emission from each nucleonic 

■ SQ —' • .. 





FIG. 7. (Color online) Neutrino emissivity via pair-breaking 
processes as a function of baryon density in units of no for 
Bie = 0 (a), Bie = 0.5 (b), and Bie = 1 (c). The full pair¬ 
breaking emissivity is shown by solid lines and consists of ^5*0 
neutron pair emission for nt/no > 0.5 and of the sum of ^So 
proton and ®P 2 neutron pair emission for larger densities. The 
separate contributions of ^5'o proton and ®P 2 neutron pairs 
are shown by dash-dotted and dashed lines, respectively. 


tify the pair-breaking neutrino emission, consider their 
neutrino emissivity, which is given hy (h = c = 1, our 
notations follow Ref. El) 


_ AG'jpim^kpn rr7 s/p 

~ 15715 


A 


s/p' 


T 


^UlpKpp rp7 ^s 


157r5 



(19) 

( 20 ) 


where the subscripts n and p refer to neutrons and pro¬ 
tons and the superscripts S and P refer to ^Sq and ^P 2 
pairing of neutrons. in Eq. (IT^ stands for the angle 
averaged value of the spin-triplet neutron gap, in which 
case it can be factored out of the integral I. (The ex¬ 
plicit form of the integral I is not needed here and can 











be found, for example, in Ref. HI). The a-coefficients 
are defined as 

anCSo) = + ^cl^vl^Xn, ( 21 ) 

«p(^‘S'o) = + -^Cpj^VppXp, (22) 

anCP2) = (23) 

where Xn/p = l + (42/ll)(m*^p/m„/p)2, c„y = 1, C„a = 
gA, Cpv = 4sin^ vv - 1, and CpA = gA, with qa 1.26 
and sin^ 9w = 0.23. 

Figure [7] displays the functions 


C. Relating the surface and crust-core boundary B 

fields 

Because only the surface Bg — 10^® G fields are ob¬ 
served in magnetars it appears to us useful to address the 
problem of relating these observed surface fields to those 
in magnetar interiors as predicted by theoretical mod¬ 
els. Equilibria of magnetized neutron stars with super¬ 
conducting cores have been constructed in Refs. [^[^. 
Both poloidal and toroidal fields, as well as their com¬ 
binations have been considered. These studies suggest a 
linear relation of the form 

Bg ~ aeBb, (25) 


Qn/p{B)[ilJl 


^n/p^Fn/Fp s/p ( ,,, , 
m^/p T ) 


which are more convenient for our analysis than the emis- 
sivities as all common factors appearing in the emissiv- 
ities m and (1701) are discarded (including the temper¬ 
ature, which is assumed to be constant throughout the 
core and the inner crust of the star). In the crust of 
the star (i.e. for densities n < 0.5no) the pair-breaking 
emission is due to the paired neutron Cooper pairs. 
This process is unaffected by the unpairing effect and is 
shown for comparison. At larger densities, in the core of 
the star, neutron and proton Cooper pair-breaking pro¬ 
cesses contribute about equally to the neutrino energy 
loss in the zero-field limit [Fig. 7 (a)]. The influence 
of the unpairing effect is seen in (b) and (c) where we 
assume constant value of the field Biq = 5 x 10^® and 
Biq = 10^® G. The constant B field removes the proton 
pair-breaking processes in the regions where B > Hc 2 
locally, because the condensate vanishes in that region. 
As a consequence the total emission rate is reduced to 
its value corresponding to the emission by the P-wave 
condensate. 

The unpairing effect will influence, apart from the 
emissivities of the magnetars, also their thermal inertia, 
because the absence of proton superconductivity will en¬ 
hance the heat capacity of the star. As a consequence 
the timescale needed for the star’s temperature to reach 
a given value will be larger than in the case of ab¬ 
sence of unpairing. In superconducting stars the main 
source of heat capacity are electrons; in magnetars non¬ 
superconducting protons will approximately double the 
heat capacity of the core of the star. Thus, we antici¬ 
pate that the proton and electron specific heats decrease 
linearly with temperature as in normal Fermi liquids, 
whereas the heat capacity of superfluid neutrons will be 
reduced by their superfluidity (exponentially in the case 
of S-wave pairing and as power-law in the case of P- 
wave pairing). The unpairing induced reduction of the 
neutrino emissivity and the increase of the specific heat 
of matter will both act to increase the cooling time-scale 
of the star. 


where Hb is the field intensity at the outer boundary of 
the core and Bg is the surface field. For purely poloidal 
field Ref. [1^ finds as = £bl^ where etR is the thickness 
of the crust, R being the radius of the star. This relation 
was derived for low fields B ^ Hd , where the role of the 
lattice of flux tubes can be neglected and these can be 
treated as isolated entities. Its validity for larger fields 
B ^ Hc 2 , more relevant to our discussion of the impairing 
effect, is not known. Nevertheless, we extracted the val¬ 
ues of £{, using our EoS shown in Sec. UTl assuming that the 
crust-core boundary is at n ~ O-Buq. We find that there 
is roughly two orders of magnitude drop in the field value 
between the crust-core boundary and the surface of the 
star. Specifically, for the 1.4 Mq star model an = 0.058 
and for the 2.67 Mq star model as = 0.021 ]28l |. The 
study of Ref. which uses a different method, sug¬ 
gests that the drop of the field from the magnetic pole to 
the base of the crust is smaller and in the limit of large 
fields is of order of unity. Clearly, further work is needed 
to establish the relation (ITS!) in the strong-field regime 
B Hc 2 - While the relation (ITS!) is highly important 
for relating the physics of the unpairing effect to the ob¬ 
servations of magnetars, our discussion and results are 
independent of the value of the coefficient as appearing 
in that relation. 


V. SUMMARY AND OUTLOOK 

We have calculated the critical field Hc 2 for un¬ 
pairing for the proton condensate, including its cou¬ 
pling to the density of the background neutron conden¬ 
sate, using Ginzburg-Landau theory in the vicinity of 
superconducting-normal phase transition. We find that 
this coupling enhances the value of the critical field by 
-30% (Eig.Hl). 

The composition of dense matter and the dependence 
of proton pairing gap on the Eermi momentum implies 
that the coherence length has a minimum as a function 
of density which translate into a maximum in the criti¬ 
cal field (see Table I). The maximum is at the crust-core 
boundary and the critical field decreases towards the cen¬ 
ter of the star. Assuming the homogeneous constant B 
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field across the core and the inner crust of the star, im¬ 
plies that magnetars with interior fields B < max Hc 2 are 
partially non-superconducting, whereas magnetars with 
B > max Hc 2 are void of proton superconductivity. 

The unpairing effect implies that the emissivity of the 
direct Urea process is only Boltzmann-suppressed due 
to neutron gap and therefore is more efficient than its 
counterpart in low-field matter, where there is an ad¬ 
ditional suppression due to proton pairing (see Figs. [5] 
and [HI which illustrate the argument in the allowed and 
forbidden kinematical domains, respectively). Unpairing 
further implies that the Cooper pair-breaking processes 
in protonic matter are absent; this reduces the local net 
pair-breaking emissivity of matter by a factor of a few 
(see Fig. III). In addition unpairing increases the specific 
heat of magnetar cores and, therefore, the thermal inertia 
of the core by a factor of two. Combined, the decrease in 
pair-breaking neutrino emissivity and the increase of the 
specific heat will increase the cooling time-scale of the 
star. This would be counterbalanced by enhancement in 
the direct Urea cooling in the strong field limit. Detailed 
cooling simulations can reveal the relative importance of 
these different factors in the cooling of magnetars, which 
we have discussed separately. 

It is not possible to state firmly whether the unpairing 
effect is operative in magnetars with observed surface 
B fields 10^® G or not, because the topology and the 
strength of the interior fields are not known accurately. If 
there is an increase of the field from the surface towards 
the center of the star (say by a factor of 10 to 15), as 
suggested by a number of studies and broadly conjectured 
in the literature, then the unpairing effect implies that 
the observed magnetars are either partially or completely 
non-superconducting. 

A separate issue, to be studied further, is the influ¬ 
ence of the strong magnetic fields on the pairing in neu¬ 
tron matter. The S-wave neutron pairing will be sup¬ 
pressed by strong magnetic field due to the Pauli param¬ 
agnetic alignment of neutron spins along the B field. The 
Chandrasekhar-Clogston limiting field for the quenching 
of S'-wave neutron superfluidity is close to the limiting 
fields compatible with gravitational equilibrium. On the 
other hand, the P-wave paired neutron fluid does not 
experience suppression in the B field [T^ . 
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Appendix: Fitting formulas 

The pairing gap for neutrons in the core and the crust 
and for protons in the core of the star were fitted by 
suitable functions (a sum of a polynomial and an expo¬ 
nential function) which depend on the Fermi-momenta of 
respective particles at zero temperature. These are given 
by 

A„(i5o) = 2.76991 - 2.17347/exp(A:p„) - 5.91497A:Fn 
-f 17.653fc|.„ - 19.1544fc|,„ -k 6.14977A:|.„, 

(A.l) 

for neutron superfluid in the crusts 

A„(3P2) = 5.97989 - 2.45018/ exp(fc|,„) - 9.76221A:Fn 
-k 6.24521A:|„ - 1.73691fc|.„ -f 0.173889A:|,„, 

(A.2) 

for the neutron superfluid in the core and using the CD 
Bonn interaction with Bruckner-Hartree-Fock spectrum, 
and 

Ap(i5o) = -302.669-f 302.982/exp(A:Fp) - 8.3717 A:fp 
- k 369.944fc|,p - 160.227fc|.p - 11.3246fc|.p, 

(A.3) 

for the superconducting protons in the core. The effective 
masses of neutron and protons were assumed to be qual 
and given by the following fit formula 

— = 1.00661 - 0.649838A:f + 0.34416fc|. 

m 

- 0.0441441A:|, (A.4) 

where kp stands for neutron or proton Fermi momentum. 
More accurate treatment would require different effective 
masses for neutrons and protons, but the corrections to 
the emissivities are expected to be small, in the range of 
a few percent. 
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